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In this paper we consider a sequential "meson emission" mechanism for charmonium decays of the 
type NNm, where is a generic charmonium state, N is a nucleon and m is a light meson. 

This decay mechanism, which may not be dominant in general, assumes that an A^iV pair is created 
during charmonium annihilation, and the light meson m is emitted from the outgoing nucleon or 
antinucleon line. A straightforward generalization of this model can incorporate intermediate A'^* 
resonances. We derive Dahtz plot event densities for the cases 'i' — rjc, J/ip, XcO, Xci and tp' and 
m = Tv'^ , fo and lu (and implicitly, any 0~^, O^""" or 1~~ final hght meson). It may be possible to 
separate the contribution of this decay mechanism to the full decay amplitude through characteristic 
event densities. For the decay subset 'I' — > ppiv^ the two model parameters are known, so we are able 
to predict absolute numerical partial widths for r(^' — > ppiv''). In the specific case J/ip — >■ ppTr" the 
predicted partial width and Afp^ event distribution are intriguingly close to experiment. We also 
consider the possibility of scalar meson and glueball searches in — >■ ppfo- If the meson emission 
contributions to — >■ NNm decays can be isolated and quantified, they can be used to estimate 
meson-nucleon strong couplings {gjvjvm}, which are typically poorly known, and are a crucial input 
in meson exchange models of the A' A'' interaction. The determination of gNNn from J/ip ^ ppiv^ and 
the (poorly known) qmnui and the anomalous "strong magnetic" coupling knnui from J/i/) — >■ ppLj 
are considered as examples. 

PACS numbers: 13.25.Gv, 13.75. Cs, 13.75. Gx, 21.30.-x 



I. INTRODUCTION 

Charmonium strong decays of the type 5" NNm, 
where \t' is a generic charmonium state, A'^ is a nucleon 
and m is a light meson, have recently attracted inter- 
est both as sources of information regarding the N* 
spectrum [H-Q and in searches for a low energy NN 
enhancement "X(1835)", which has been reported in 
J/ip -J> jpp [i] and J/'0 jn+TT^if 0, H, but thus 
far not in ^ ^ ppm. These decays are also of interest 
because their partial widths can be used to estimate the 
pp associated charmonium production cross sec- 

tions at PANDA d, [l^. As we shall show here, they 
may also provide information on NNm meson-nucleon 
coupling constants, which could be used to identify un- 
usual resonances such as molecule or glueball candidates. 

Specific ^' NNm reactions that have recently been 
studied experimentally include J/^p — pp7r° [ill 

PPV 



and pp-q' |4|, and ppujM; ip' pvir' 
pnir^ -t- h.c. 
limit) pp(j) [ill . 

These decays may prove to be complicated processes in 
which several decay mechanisms contribute significantly. 
For this reason it will be useful to have predictions for 
\1/ — > NNm Dalitz plot (DP) event densities assuming 
various decay mechanisms; this paper provides these re- 




ppp [llj and ppoj [ill, |l^ , and (upper 
and XcJ ppn^ and pprj [iT 



suits for one such mechanism. In particular we derive 
the DP event densities that follow from sequential meson 
emission, in which the charmonium state (generically ^E*) 
decays to an intermediate NN state, which radiates the 
light meson from the N or N line, — )■ NN NNm. 
The two Feynman diagrams assumed in this model are 
shown in FiglTJ 

We emphasize that the actual relative importance of 
this and other ^ — ^ NNm decay mechanisms is unclear 
at present, and may depend strongly on the charmonium 
state ^' and the light meson m; one purpose of this paper 
is to determine the rates predicted by this meson emission 
decay model in isolation for comparison with experiment, 
so that the importance of this decay mechanism can be 
estimated. 

The predictions of this NN — ^ NNm decay 

model can be given in some cases with no free parameters, 
since the strengths of the a priori unknown couplings 
'ilNN and NNm can be estimated from other processes. 
Here we will give absolute predictions for the set of par- 
tial widths {r(^' -> ppTT°)}; we use the known partial 
widths {r(^' —5- pp)} to estimate the {^'A^TV} couphngs. 
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FIG. 1: Feynman diagrams of the meson emission model. 
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and the final NNtt coupling is of course well known. 

Provided that the contribution of the 4* NN — > 
NNm decay mechanism can be isolated and quantified 
experimentally, this information can be used to estimate 
meson-nucleon strong coupling constants; these are gen- 
erally poorly known, and play an important role in nu- 
clear physics as input parameters in meson exchange 
models of the NN force [M^^]- 

We will also discuss the determination of the (well 
known) ppn^ and the (poorly known) ppuj couplings from 
the decays J/ip — > PpT^^ , ^' ~> ppT^^ and J/ip pp^^ as 
examples. This provides a third motivation for the study 
of ^ — ?> NNm decays; they may prove useful for estimat- 
ing NNm coupling constants, in addition to their rele- 
vance to N* spectroscopy fT-o| and low-mass pp dynam- 
ics [^-Q. Another motivation for studying ^E* — )> NNrn 
is the possibility of observing light scalars, including the 
"a", the 980 MeV states and the scalar glueball, in the 
decays 4* — >■ ppfo (and Oq). 



II. FORMULAS 



The areas {^_d} of these Dalitz plots are useful for es- 
timating r(^ ppm) partial widths i9|. Although Ad 
can be evaluated in closed form for \1/ — >■ ppm with gen- 
eral mass ratios r and i?, the resulting expression is quite 
lengthy, so when required we will simply evaluate each 
Ad numerically. 

In deriving the DP event densities we have usually as- 
sumed that the '^pp coupling is a constant times the 
simplest relevant Dirac matrix for the given ^ quantum 
numbers; for example, for the J/^ we use a pure vector 
J/ipPP vertex, —igj/^ppjfj,- The order of the hadron la- 
bels in gabc is meant to reflect the fact that the numerical 
value of this coupling constant is taken from an a — > 6c 
transition, here J/ip ^ pp. This could be a significant 
concern if form factor effects are large. 

We proceed similarly for the light mesons tt" and /o; 
for the pion we use a pure pseudoscalar NNn coupling, 
with vertex (7Arjv7r75, and —igNNfal for the NN/q vertex. 
Since light vector mesons (generically represented by the 
w) have two interesting strong couplings, Dirac (vector) 
and Pauli (anomalous magnetic) , for this special case we 
assume a vertex with two interactions. 



Here we will usually specialize to charmonium decays 
to a pp pair and a neutral meson, vj/ ppm^; these 
decays are reasonably well studied, and enjoy the simpli- 
fication of equal baryon and antibaryon masses. Our re- 
sults employ conventions for kinematic variables, meson- 
baryon couplings and masses that were used in Ref. "24^ . 
In particular, M is the mass of the initial charmonium 
state, mp is the proton mass, rn™ is the mass of meson 
(subscript) m, and dimensionless mass ratios R = M/mp 
and r = mm/mp are defined relative to the proton mass. 
(Hence the numerical values of R and r depend on the 
decay process.) We use scaled dimensionless variables 
X = Mp^/m^ — 1 and y = Mp^/m^ — 1 and their in- 
verses u = 1/x and v = 1/y io describe Dalitz plots; 
these greatly simplify our results. The DP event densi- 
ties we present here are formally partial width densities 
in X and y, which are related to the more familiar differ- 
ential partial widths by a trivial overall constant. 



p(Lj) 

A' 



igNNu ( 7ai + i{i^NNLo/'^mp)a^uqv 



(5) 



Ref. [23| assumed a similar J/ipPP vertex; see Ref. [2j for 
additional details regarding the couplings assumed here. 
We generally abbreviate these coupling constants as g^ = 
9^1 NN gm = gNNm] rationalized squared couplings 
= g|,/47r and «,„ = g^/47r are also used. 
For the special case of ppir^ final states, these event 
densities can be obtained by applying crossing relations 
to our previously published results for the unpolarized 
differential cross sections for the 2 — ;> 2 processes pp — J> 



TT^Vl/ [2J|; the other (/o and w) cases have not been con- 
sidered previously. The results for all cases considered 
here are given below. 



(PT{ric -J> pp7r°) 
dxdy 



{u 



d'^TNj ^ ppm) . d'^T{'^ ^ ppm) 

^ m^ T^TT; T^TT^ . (1) 



dxdy 



dM^^dM?„ 



Before we give our results for these event densities, it 
is useful to recall some general properties of a ^ — t- ppm 
Dalitz plot. The boundary in the dimensionless variables 
(x, y) is specified by the curves 

r2i?2 + (^2 + _ 2)x -x^ ± F,nF^ 

^ 2{iT^) 

where Fm — F{r, x) and F^ = F{R, x), with 

F{a, x) = (a2(a2 _ 4) _ 2a^x + x^/^ (3) 
The range of values of x (and y) in the physical region is 
r{r + 2)<x<R{R-2). (4) 



--r2i?2 (6) 



d'^r{ J/ijj ppT:°) _ mp j {u + w)2 

-2{u + v){u + v + l)r'^ + 2uvr'^ - (u^ + v"^) ■ 
r2i?2l (7) 



TTlp I . ,9 



dxdy 

f^+4r2-r2i?2^ 
\uv ) 



(8) 
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(fT{xci pp'K°)_ _ rup I (u + v) 



dxdy 



(9) 



d^T{J/^^ppfo) 



(10) 



uv 



+8{u + v){u + V + 1) -2^u{u+l)+v{v + l) 

+6uvj + 4(m + vf R^ + 2uv r"^ - {u^ + v"^) ■ 

dxdy 87ri?3 \ uv 

-16(m + v){u + v + l)+4{u + vf{r^ + R^) 

-{u-vYt'^rA 



(11) 



(12) 



<P^{Xci PPfo) 
dxdy 



{u + v) 



dTrR^ 



uv 



{u + v + l) + r^+(^^^--^-8{u + v)- 
\ Zuv 

{u + v + 1)) i?^ + (2{u^ + v^) + 8uv + u + v^ r^R^ 



2d4 „,„,4d2 {u'^+V^) 



-2{u + vy R^ -uvr^R^ - 



d'^V[ijc ppuj) m 
- = a, 



(13) 



(m + 1;)^ 



dxdy 47ri?3 I [ uv 

2{u + v)iu + v + l)R^ - (w^ + v^) r^R^ + 2uv R'^ 

2{u + vf 



uv 



+ {^{u^ + v^) + 2uv^ r^R^ 



{u + v){u + V — 1) {u + v) 



2uv 



+ 



8uv 



■r' + ^R^ 



+ 



^^-{u^+v^)) r^R^-^^^^r'R^ 



--T^R'' 
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(14) 



d'^T{J ppuj) m 
— = J/v 



+ w^) 



dxdy " 6^i?3 

-4(u + v){u + V + 1)- 2{u{u + 1) + v{v + 1)) • 

(r^ + R^) + 2uw - {u^ + - 4uw) r^ii^ + 2uv R^ 

(u -\- v)^ 

+Ku, \- 6{u + v){u + V + 1) — 6uvr'^ 



uv 



-(3{u^ + v'^) -8uv^ r^R^ 



{u + v) 
2uv 



{u + vf 
8uv 



- 2{u + v){u + v + l)^r^ + ^R^ 



-^(u{u + 1) + v{v + 1) - 6uv^ - (u - vf r^R^ 



+ -uvr^ - -(u^ + u^ -Auv)r'^R^ 
4 8 



(15) 



{u + vf 
uv 



rf^r(xco -» ppuj) 

dxdy 

+8{u + v){u + 1; + 1) + 4(m + vf - 2{u{u + 1) 
-\-v(v + 1) + Quv^ R^ - {u^ + v^) r^R"^ + 2uv R* 

- 12{u + v){u + v +^)r'^ + 3{u + vf r^R^ 

{u + v){u + v + l) f {u^ + 6uv + v"^) 
2uv \ 8uv 

+^^u + v){u + v + \)\ r2 + ii?2^l!i±^r^ 



u^ + Auv + v'^ + 



{u + v) 



r^R' 



(u — v)^ 



+ ^r'R' 



rf^r(xci ^ ppuj) 

dxdy 



r^R^ 
(16) 



rur, 



{u + vf 



uv 



+ 



(^^^^^p- + 4{u + v){u + v + l)^ R' 



+ (2{u -vf -u-v^ r'^R^ - (w^ + v"^ + Quv 



+u + v)R'^ + uv r^R"^ + (2uv - 
+uvR'^ 



{u + vf 2 ("^ + - Auv) 

UV 2uv 



R^ - (^{u^ +v'^)+u + v^ r^R^ - UV r'^R^ 

'{l + 2{u + v) + 2{u + vf) 



-\{u'+v^)t^R' 



Auv 



4 



{u + V + {u - v)'^) 2 /3 (u + v) 



8uv 

.2 I „,2 



1 



+ 2{u^ + v- uv) + 



2 4mu 
(u + v) {u^ + v"^ + Auv) 
2 l&uv 



r^R^ + -R'^ + -iu + v + 2{u^ + + Quv) ] r^R^ 
4 8V ' ' 

2^ ' 8 



1 

16 



[u^ + + Auv) r^R"^ 



(17) 



The symmetry of these event densities under {x, y) and 
hence {u, v) interchange is a consequence of C-parity in- 
variance. There are singularities in these events densities 
along the hues x = G {u = oo) and y = Q [v = oo), corre- 



sponding to M 



2 

pm 



ml and M?„ 



These are due 



to the poles of the p and p propagators in the Feynman 
diagrams for the decay process \1/ — )■ pp — )■ ppm (FiglT|) 
and he outside the physical regions of the Dalitz plots. 
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FIG. 2: The J/xj) ppTv'^ DP event density predicted by the 
meson emission decay mechanism J/tp pp ^ Pp!^'^, Eq.lO. 
Contours of equal density are shown. 



III. APPLICATIONS 

A. The J/ip — > pptt" partial width 

As a first application we will evaluate the partial width 
for J/ijj — ppn'^ assuming this meson emission decay 
mechanism. The PDG [26] quotes a branching fraction 
of 



B{J/iP pPtt^) = 1.09 ± 0.09 • 10 



-3 



(18) 



which is the average of early measurements by Mark-I 
ll], BASF [III and Mark-II There are also more 

recent experimental results on this decay from BES-II [l| . 
Using the current PDG value for the J/4> total width of 
93.2 ± 2.1 keV, this branching fraction corresponds to a 
partial width of 



F( J/V- ppTT^) = 102 ± 9 eV. 



(19) 



To evaluate this partial width in the meson emission 
model we simply integrate the theoretical event density 
([7]) over the Dalitz plot. This event density is given by 



d^r{J/Tp -> pptt") 
dxdy 



aj/^a^mpp{x,y) (20) 



where the dimensionless density function p is 
1 f {u + vY 



pix,y) 



127ri?3 



2{u + v){u + V + l)r^ 



+ 2uvr'^ - {u'^ + v^) R' 



(21) 



This (parameter-free) relative event density, scaled to the 
maximum value in the physical region, is shown in Figl^l 
Integrating (pn|) over the Dalitz plot gives 



T{J/i> ppTT°) = aji^ a.„mp ■ (p) ■ Ao/mp . (22) 

where (p) is the mean value of p{x, y) in the Dalitz plot, 
which has (physical) area A^; 



DP 



pdxdy = (p) ■ Ad I nip . 



(23) 



We evaluate (p) and Ab numerically, assuming phys- 
ical hadron masses; we use FDG masses rounded to 
0.1 MeV; m^o = 0.1350 GeV, vip = 0.9383 GeV and 
TO,7/.0 = 3.0969 GeV, which leads to (p) = 3.070 • lO^^ 
and Au — 9.265 GeV'', and a partial width of 



F( J/V' ^ ppT:°) = 34.44 • a 



MeV. 



(24) 



To complete this estimate we require numerical values 
for the NNtt and J/ipNN coupling constants. For NNtt 
there is general agreement from meson exchange mod- 
els of NN scattering that gNNn ~ 13 (see for example 
[T^ - [2l| ): we accordingly set gNNir = 13.0. The value of 
the J/ipNN coupling constant (here gj/ippp) can be es- 
timated from the measured partial width to pp, which is 
(again using PDG numbers) T[J/^ — )• pp) = 202 ± 8 eV. 
Equating this to the theoretical decay rate 

V{J/^~^PP)^ aj/^/3p(l + 2/i?2)M/3 (25) 

gives a value of qjhhf,p — 1-62 • 10^'^, as was quoted previ- 
ously in Ref. 2J| . Using these couplings in Eq. (|24| gives 
our meson emission model prediction 



T{J/i!^ pp7r°) = 97 eV. 



(26) 



This is consistent with the experimental value of 102 ± 
9 eV quoted in Eg. p^ . This excellent agreement is 
somewhat fortuitous, since this version of the model 



does not include the N* contributions evident in the 
J/ip ppn° Dahtz plot [l| (see also FigE]). 

We note in passing that the charged-pion cases J/ip — J> 
pnT:~ and npir^ should have branching fractions close 
to twice B{J/iIj ppn'^), reflecting an isospin factor of 
two. Experimentally this is indeed the case; the ratio of 
the PDG J/iJj branching fractions to pn'K~ and ppTr" is 
B{J/ij ^ pnTT-)/B{J/i; ^ ppTT^) = 1.94 ±0.18. 

B. Projected J/ip ppyr" event densities 

Projections of DP event densities onto the invari- 
ant mass of one pair of particles are useful in searches 
for intermediate resonances. For J/ip — > ppTr'^, the 
BES Collaboration has published acceptance-corrected 
event densities in Mp^^ and Mp^^ invariant mass (Fig. 6 
of Ref.jH), which show clear evidence for N* reso- 
nances. Here we will generate the corresponding theo- 
retical Mp^-projected event distributions in the meson 
emission model for comparison with experiment. Al- 
though N* resonances are not incorporated in our calcu- 
lation, this comparison will test the relative importance 
of the meson emission decay mechanism in this decay, and 
establish whether the model predicts a non-iV* "back- 
ground" invariant mass distribution that is similar to the 
data in form and magnitude. 

The full two-dimensional DP event density (PT/dxdy 
predicted by the meson emission model is given by Eq. (jT]) . 
Projecting this onto Mp^^ is straightforward; first one in- 
tegrates over y = Mp^/nip — 1 between the DP bound- 
aries y±{x) of Eq.®, which gives dT/dx. Converting this 
into a distribution in Mp^r introduces a Jacobean, which 
is specified by the definition x — Mp^/rri^ — 1. This gives 

dT{J/iP ppiT^) _ 2Mp^ d^T{J/ij pp7r°) 

dMp^ Tfip Jy_(x) dxdy 

(27) 

We have evaluated this distribution numerically, given 
the d'^T{J/ip ppn")/dxdy of Eq.© and the masses 
and couplings a,//^ and a-rr used in Sec. III. A. The result 
is shown in Figl21 together with an experimental distri- 
bution provided by BES (reported in Ref.[l|), using a 
common scale. The data is the combined acceptance- 
corrected Mp^ and Mp^ distribution, scaled to give their 
reported B{J/il} ppn^) — 1.33 • 10^^ rather than the 
PDG value of 1.09 • 10"^. 

Clearly there is a close resemblance between the me- 
son emission model prediction for the J /ijj — pp7r° event 
distribution in Afp^ and the observed BES distribution, 
both in form and magnitude. This suggests that a study 
of this reaction assuming this model for the experimen- 
tal "background" combined with N* resonance contribu- 
tions would be an interesting exercise. Although BES [ij] 
recently reported a similar study, they introduced an ad 
hoc St^n(= Mp^)-dependent form factor that suppressed 
this "background" meson emission amplitude relative to 
N* contributions. The similarity to experiment evident 
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FIG. 3: J/V' — >■ pp7r° experimental (BES) and theoretical (me- 
son emission model, Fig[2} Dalitz plot (DP) event densities, 
projected onto Mp^. This is not a fit; see text for discussion. 

in Fig|3] suggests that this mechanism merits additional 
consideration. 

C. Other — > ppTr" partial widths 

Since the two meson emission model parameters 5*pp 
and gNNm are both known for several \1/ — >■ ppn^ decays, 
we are able to give absolute predictions for these partial 
widths. (We previously used F(^ — pp) to estimate 5*pp 
[131; here we again use these values, and set gNNyr = 13.) 
These ^ ppn^ partial widths are given in Table HI 
together with some intermediate theoretical quantities 
and the experimental widths. (These experimental values 
are the PDG total widths times branching fractions, with 
errors added in quadrature.) 



* 


10'' ■ (7*pp 


10' ■ ip) 


Ad [GeV*] 


- n 
pptt" 


- n 




19.0 


0.530 


6.862 


1.7 keV 




J/i, 


1.62 


3.070 


9.265 


97 eV 


102 ± 9 eV 


Xo 


5.42 


3.691 


18.605 


2.6 keV 


6.0 ± 1.3 keV 


Xi 


1.03 


0.554 


22.351 


17 eV 


103 ± 43 eV 


^' 


0.97 


2.010 


30.501 


75 eV 


41 ± 5 eV 



TABLE I: Comparison of theory (meson emission model) and 
experiment for r('I' — ppyr") (see text). 



These rates were derived using Eq.([25]), with the ap- 
propriate density function p chosen from the set Eqs.(IB]- 
ini). In addition to the rates. Table U also gives the cou- 
pling constants assumed, the average of the density func- 
tion p over the Dalitz plot, and the DP area Ab in phys- 
ical units. 

It is clear from the table that the wide variation in the 
absolute scale of partial widths observed experimentally 
is approximately reproduced by the model, at least at 
a "factor-of-two" level of accuracy. This suggests that 
the meson emission decay mechanism may indeed be an 
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important component of the decay amphtude in all these 
decays; a more definitive test would involve a direct com- 
parison of the DP event densities or their two-body pro- 
jections, as in FiglJl 

The Xci case appears to be an exception to this approx- 
imate agreement, however in view of the large experimen- 
tal error it is not clear if this is a real discrepancy; theory 
and experiment only differ by 2a. 

Although the single experimentally unobserved decay 
Vc ~^ pp'n:^ is predicted by the meson emission model 
to have a relatively large partial width of 1.7 keV, it 
is actually considerably suppressed by the presence of 
an on-diagonal node in the DP event distribution. An 
experimental study of rjc — > ppi:^ would accordingly be 
very interesting, since the contributions of other decay 
mechanisms may be easier to identify in the region of the 
DP where the meson emission model gives a zero or weak 
contribution. 



D. qnnt! from B{J/tp ppn'^) / B{J/4^ — )► pp) 



tp' case has a much larger ppir^ DP area than the 
J/V', which is partially compensated by a smaller mean 
event density (p). (These quantities are given in Ta- 
ble m) Using M = 3.6861 GeV and the PDG branch- 
ing fractions B{ip' -J> pp-K^) = (1.33 ± 0.17) • 10"^ and 
B{ip' pp) = (2.75 ± 0.12) • 10"^, we find the V'-based 
gNNiT estimate 



ONNtt 



= 9.9 ±0.7 . 



(30) 



'0'— >pp7rO 



This is similar to but somewhat smaller than the estimate 
obtained above from Jf-ip decays (|29p. and may give an 
indication of the accuracy of this approach for estimating 
NNm coupling constants. 

Of course the other relations for am analogous to (pS)) 
will only be useful if the meson emission decay mech- 
anism is dominant in those decays as well. Otherwise 
the contribution of this mechanism to the decay must be 
identified and quantified, for example through a detailed 
study of the DP event density. 



Previously we noted that 4' — ppm decays can be used 
to estimate NNm couplings, provided that the contribu- 
tion of the meson emission decay mechanism to the de- 
cay amplitude can be quantified experimentally. In the 
following we will use the decay J/ip — ppn^ as an il- 
lustration of this approach, since the agreement between 
the experimental and theoretical partial widths suggests 
that domination of this decay by meson emission is a 
reasonable first approximation. 

Since the a priori unknown coupling aj/^ cancels 
in the theoretical branching fraction ratio B{J/ip 
ppTT^)/B{J/ip pp), we can use it to estimate gNN-n 
directly. The meson emission model decay width for 
T{J/ip ppTr°) and the two-body decay width (1^ 
imply the following relation between this ratio and the 
coupling = g^^^/47r: 



^ n-4/R2U/2 (R + VR) B{J/^P^ppn°) 
^ ^ Hp)An/m^^' B{Jl^^pp) 



(28) 



Substitution of the experimental PDG numbers 
B{J/iP ppiT^) = (1.09 ± 0.09) ■ 10-3 and 

bIj/^P pp) = (2.17 ± 0.07) • 10-3 for these branching 
fractions leads to the estimate 



QNNtt 



= 13.3 ±0.6 



(29) 



which is consistent with NN meson exchange model val- 
ues. 

We expect to find approximately equal gNN-n estimates 
from other ^ pp and ppm decay pairs if the me- 
son emission decay mechanism ^ ^ pp ^ ppir^ is in- 
deed dominant. A second state ^' that can be use to 
estimate gNNn is the i/)'(3686). Since the ip' has the 
same quantum numbers as the J/ip, Eq. (j28p is again 
appropriate for our coupling constant estimate. This 



E. Scalar mesons in *!/ 



NNm 



The long-standing interest in the light scalars makes 
the possibility of studying them using these decays an 
attractive prospect. This motivated our inclusion of de- 
cay formulas for the processes ^' — >■ ppfo in our set of 
theoretical DP event densities. 

Here we will give meson emission model predictions 
for the branching fractions B{'i' PPfo), where 5* = 
ric J/iPtXcO,Xc1 and ip', for a light "sigma" meson with 
TO/o = 0.5 GeV. To evaluate these partial widths we 
proceed as in Sec. III. C, and integrate the appropriate 
decay width formula from the set Eqs. (|10m3p over the 
Dalitz plot. We again use the "ifpp coupling constants of 
Sec. III. C, as given in Table HI The total widths used to 
convert the calculated partial widths to branching frac- 
tions are the current PDG values, r(»7c) = 27.4 MeV, 
r(J/V') = 93.2 keV, r(xco) = 10.4 MeV, r(xci) = 
0.86 MeV and ^{ip') = 309 keV. Our results are given 
in Table Hll Since there is no general agreement re- 
garding an NN fo{500) coupling constant, in the ta- 
ble we first give the predicted branching fraction rela- 
tive to B{J/^p — >• pp/o(500)) = Bq, which is numeri- 
cally 0.338 ■ 10"'' ■ g'jqpff^- (The unknown gNNfo can- 
cels in these ratios.) The next column gives absolute 
i3('I' — > ppf o{500)) branching fractions for a rather arbi- 
trarily chosen gNNfo = 10- Finally, the table quotes ex- 
perimental branching fractions for the related processes 
^' ppTr+Tf- for comparison. 

The relative theoretical branching fractions (Table [TTl 
col. 2) suggest that the best channel for identifying a light 
scalar in \1/ — ppfo is J/ip ^ ppfo itself (assuming that 
the meson emission model is a reasonable guide) . Given a 
somewhat larger event sample, ip' — >■ ppfo should be com- 
petitive with J/ip, and has the advantage of more phase 
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space, so the scalars near 1 GeV and the /o(1500) could 
also be investigated, rj^ — >■ ppfo has a comparable the- 
oretical branching fraction, but the difhculty of produc- 
ing the rjc makes this a less attractive channel. Finally, 
the XcJ states are predicted to have much smaller ppfo 
branching fractions than J/ip — >■ ppfo, and accordingly 
are less attractive experimentally if this decay model is 
accurate. 







10'^ • ]if^y-^3=w) 

ppfo 


jq3 . ^e.pt. 




0.40 


1.4 


< 12 (90% c.l.) 


J/V 


= 1 


3.4 


6.0 ±0.5 


Xo 


0.045 


0.15 


2.1 ±0.7 


Xi 


0.016 


0.054 


0.50 ±0.19 




0.21 


0.72 


0.60 ± 0.04 



TABLE II: Theoretical (meson emission model) branching 
fractions for light scalar meson production. The numeri- 
cal columns are 1) The ratio _B(* -> pp/o(500))/B( J/^ 
pp/o(500)); 2) 10^ ■ B(* ^ pp/o(500)) for gNNjo = 10; 3) 
lO^ . B'=^P* (^» ppTT+Tv'), for comparison with 2). (See 
text.) 

We have also estimated the effect of an /o(500) width 
on these results. Imposing a Breit-Wigner /o mass profile 
with T-fg — 0.5 GeV, truncated at 2171^^, decreases all the 
absolute theoretical pp/o (500) branching fractions 
in Table |TT] (col. 3) by « 10%. The relative theoretical 
branching fractions (col. 2) are even less sensitive to the 
/o(500) width, and become 0.41,0.048,0.017 and 0.21. 

A light scalar fo meson would presumably decay 
strongly and perhaps dominantly to tttt, so decays of 
the type ppmr are of special interest, notably 

J/ip — >■ ppTTTT (in view of our large theoretical B{J /ip — > 
pp/o(500))). The charged case JN) — >■ pp7r+7r~ has been 
studied by Mark-I 27], DESY [3Q| and Mark-II [2|. Al- 
though no light scalar mesons have yet been identified 
in this decay, it is suggestive that J/ip — > pp'!T~^Tr~ is 
the largest known exclusive J/^ — >■ ppX mode, with a 
branching fraction of B{J/ip — > ppn^Tr^) = (6.0 ± 0.5) • 
10-3. 

In addition to the ppfo intermediate state of interest 
here, this decay may also receive important contributions 
from NN*, N*N* and AA, as well as other two-baryon 
and NNm states; this may complicate the comparison 
with experiment considerably. Ref.[2^, which has the 
largest event sample, finds a large but not dominant AA 
contribution, B{J/tp A++A--) = (1.10±0.29) • lO^^, 
and gives a rather tight upper limit of ss 5% on the con- 
tributing subprocess Jf-ip pp(f\ B{J/ip PPP^) < 
3.1 • 10-4 (90% C.L). As Ref.[2g| shows in their Fig.31 
that the tt'^tt~ invariant mass distribution from non-AA 
events is a broad sigma-like distribution, there may well 
be a large J/^ PPf oi'^ 500) ppir'^TT^ contribution 
to this decay, with a branching fraction comparable to 
the theoretical 3.4 ■ 10"'^ predicted for gNNfo = 10 (see 
Table |TT| . It will be very interesting to investigate this 
possible light fo contribution in a future experimental 



study, as well as to search for the /o(980) and ao(980) 
scalar states and the scalar glueball candidate /o(1500) 
in (higher-mass) charmonium decays, notably of the V''- 



F. Decays to NNuj and NNV 

Charmonium decays to NNto are especially interest- 
ing, since the w plays a crucial role in meson-exchange 
models of the NN force, as the dominant origin of the 
short-ranged "hard core repulsion" , through t-channel uj 
exchange. Conceptual problems with this w-exchange 
mechanism include 1) the very small NN separation im- 
plied by this mechanism {Rnn ~ l/my ~ 0.3 fm), at 
which quark-gluon dynamics may be a more appropriate 
description of the interaction, and 2) the prediction of 
a corresponding short-ranged NN attraction and deeply 
bound NN states, which are not observed. (See for ex- 
ample Refs.|3ll-[33|. and references cited therein.) 

There are two strong coupling constants in the NNut 
vertex, as summarized by Eq.(I5]), the overall strength 
gNNuj of the Dirac (7^) coupling, and the relative strong 
magnetic Pauli coupling knnui (here abbreviated g^^ and 
K^, with a^j = g'^/Air). Unfortunately the NNV cou- 
plings in the meson exchange models are not a priori well 
established experimentally, and are therefore usually fit- 
ted directly to NN scattering data. These NN scattering 
studies are thus fits to the data rather than predictions 
that test the assumed vector-meson-exchange scattering 
mechanism. These NN fits typically find go, « 10-15 for 
the Dirac NNuj coupling, whereas the NNuj Pauli cou- 
pling has remained poorly determined; examples of NNoj 
parameter sets from the NN scattering literature in- 
clude (g^,K^) = (12.2,-0.12) (Paris), (12.5, ±0.66) (Ni- 
jmegen), and (15.9,0) (CD-Bonn) (these are cited in 
Ref . [34| ) . Independent estimates of the NNuj coup ling 
from experiment have been reported by Sato and Lee [35| 
(from pion photoproduction) and by Mergell, Meissner 
and Drechsel 36] (from nucleon EM form factors). Sato 
et al. assumed k^^ = 0, and quoted the range g^i = 7-10.5 
for experimentally favored values of the Dirac coupling. 
Mergell et al. found a small k„ but a much larger Dirac 
coupling, igu,,K^) = (20.86 ± 0.25, -0.16 ± 0.01). The- 
oretical calculations include a QCD sum rule result of 
Zhu [13, who finds {g^, n^) = (18 ± 8, 0.8 ± 0.4), and a 
recent '^Pq quark model calculation of NNm couplings 
[3^ 1 which found the analytic result — —3/2. 

Charmonium decays to ppuj final states (and NNV 
more generally) may allow independent estimates of these 
NNuj and NNV couplings, again provided that they are 
dominated by the meson emission decay mechanism, or 
at least that this contribution to the decay amplitude can 
be isolated and quantified. In the following discussion we 
will consider the decay J/tp ppuj as an example. 

Results from experimental studies of the decay J/ ■0 
ppuj have been published by Mark-I [23| , Mark-II [2^ and 
most recently BES-II 5]. The PDG value of the J/tp ^ 
ppuj branching fraction, estimated from these results, is 
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B{J/'ip — J> ppuj) = (1.10 ± 0.15) • 10~^, which combined 
with the PDG J/ip total width gives an experimental 
partial width of 

r( J/V- ^ ppuj) = 103 ± 14 eV. (31) 

The fact that this is approximately equal to the ppn'^ 
partial width (IT9)) despite the much smaller phase space 
suggests a robust NNuj coupling. 

On evaluating this theoretical decay rate by in- 
tegrating Ea.([T5|) numerically with physical masses, 
both NNuj couplings free, and (as used previously) 
9j/i(>pP = 1-62 • 10"^, we find 

r( J/V' ^ ppuj) = ■ (2.468 - 1.101 + 0.886 K^j eV. 

(32) 

The single number V{J /ip — ;> ppw) alone does not suf- 
fice to determine the NNuj strong couplings since there 
are two free parameters, and k^i- If we set k^^ = 0, 
following CD-Bonn [l^ and the Sato-Lee photoproduc- 
tion study fss'], the measured partial width (|31l) and the 
theoretical decay rate ([5^ imply 5^ = 23 ± 3 (provided 
that meson emission dominates this decay). This guj is 
rather larger than these references prefer for ; it is con- 
sistent however with the EM form factor fit of Mergell et 
al. [3^ and the range 18 ± 8 reported by Zhu ^ from 
QCD sum rules. If we instead assume the '^Pq quark 
model value — —3/2 for the Pauli term we find 
g^ = 14.6±2.0, which is consistent with the values quoted 
by NN scattering studies, and is somewhat larger than 
the photoproduction value. 

It is possible to estimate the two parameters g^^ and k^j 
independently through a more detailed comparison be- 



tween J/-0 — ppuj data and the theoretical DP event den- 
sity, Eq. p3| . This theoretical event density is strongly 
dependent on the Pauli coefficient Kj^; with = the 
event density at lower Mpp is strongly suppressed (see 
FigH. 

In contrast, for moderately large \ku\, such as the 
quark model value — —3/2, the theoretical DP 
event density is closer to uniform. This is illustrated in 
FigO which shows this event density along the diagonal 
^'^puj = ^puj (relative to the maximum value on diago- 
nal) for various < 0. If the meson emission model does 
give a good description of this decay, evidently it may be 
possible to determine Koj by comparing the J/?A — ?> ppuj 
DP event density on diagonal to the prediction in FiglS] 
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FIG. 5: Theoretical J/ili ^ ppu DP event density along the 
diagonal M^^^ = Mp^, showing the strong k^j dependence. 
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FIG. 4: Theoretical relative J/ip — > ppu) DP event density for 
= 0, from Eq. (|15p . Note the suppression near pp threshold 
(upper right). 



G. Other ^ — > ppV decays 

Other decays to ppV final states that are closely related 
to J/ip ppi^ include J/ip ppp'^ and J /if) Pp4>- 
In the meson emission model these are both described 
by the decay formula ([T5|) . albeit with different vector 
meson masses and NNV couplings. Given the rounded 
PDG masses rupo = 0.7755 GeV and = 1.0195 GeV, 
we predict the numerical decay widths 

T{J/^ -J. pPp°) = ap ■ (2.6I4 - 1.155 Kp + 0.930 nfj eV 

(33) 

and 

F( J/V' pp(j)) = ■ (0.184 - 0.109 + 0.087 nfj eV. 

(34) 

The /9° case is especially interesting due to the range of 
values reported for discussed by Brown and Mach- 

leidt [11]. Although vector dominance predicts Hp = 3.7 
"weak p" , and some data has been interpreted as sup- 
porting this, fits to TTTT — > NN and S-D mixing in NN 
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scattering prefer a larger value "strong p" ; the Bonn po- 
tential model for example uses {ap,Kf^ = (0.84,6.1) tlj. 
A QCD sum rule calculation by Zhu [39| finds {gp, Kp) = 
(2.5±0.2, 8.0±2.0), comparable to the fitted Bonn values. 
In contrast, the valence quark model with a ■^Pq NNp 
coupling predicts a much smaller Kp — ~Kui — +3/2 [34| . 

Using Eg 1331 we can give meson emission model predic- 
tions for r{J/ip ppp'^) that follow from these various 
{gp,Kp) parameters. The Bonn parameters cited above 
give T{J/il; -J> ppp°) = 25 cV and B{J/tp ppp°) = 
2.7 • 10~^; this is essentially equal to the current ex- 
perimental upper limit [H, S of 3.1 ■ 10"'' (90% c.L), 
which is a Mark-II result dating from 1984. The Zhu 
QCD sum rule central values for (g^, Kp) give essentially 
identical results. In contrast the valence quark model 
with '^Pq couplings gp = .9^/3 and Kp — +3/2 (and us- 
ing gp = 14.6 from the J/^ — )■ ppu discussion above) 
gives a much lower T{J /ijj — >■ ppp^) = 5.6 eV and hence 
B{J/^ — J> ppp°) = 6.0 • 10~^, which is about a factor of 
5 below the current experimental limit. The proximity 
of the Bonn and QCD sum rule parameter predictions 
for B{J/ip PPP^) to the current limit suggests that an 
experimental study with significantly improved sensitiv- 
ity could make a useful contribution to establishing NNp 
couplings. 

The decay J /ip pp4> in contrast has been observed, 
and has an experimental (PDG) branching fraction of 
Qexpt.f^ jl^ pp^-^ _ j-^ g _|_ g-j . corresponding 

to V^'P^-iJ/ip ppcj)) = 4.2 ± 1.4 eV. Unfortunately 
in this case we do not have a theoretical estimate for ei- 
ther or since the (valence level, leading-order) ^Pq 
model predicts no NN(p coupling. Clearly it would be 
very interesting to obtain experimental values for these 
couplings, since little is known about the properties of 
Zweig-suppressed vertices. Again, if the meson emission 
model gives a good description of this decay, a compari- 
son of Eq. p^ to the experimental J/ip —?■ ppcj) DP event 
distribution should allow an experimental determination 
of the NNip couplings. 

Finally, we note in passing that -0' decays to NNV 
are apparently not in agreement with the meson emis- 
sion model; proceeding as above, we would predict a 
branching fraction of B{ip' — >■ ppui) = 9.4 • 10"'', whereas 
the PDG experimental value is an order of magnitude 
smaller, B'^'^p*- {ip' ppu) = (6.9 ± 2.1) • 10"^ Possi- 
ble explanations for this discrepancy, including form fac- 
tors and (destructive interference with) additional decay 
mechanisms, are discussed in the next section. Since the 
total ip' — » ppuj data sample reported by CLEO [ll| and 
BES 1131 comprises only about 35 events, it is not yet 
possible to establish the reason for this large discrepancy 
between experiment and the meson emission model. This 
would ideally involve a comparison between the predicted 
and observed Dalitz plot event densities. Hopefully this 
comparison will be possible using the large '0' data set 
being accumulated at BES-III. 



IV. SUMMARY, CONCLUSIONS, AND 
FUTURE DEVELOPMENTS 

In this paper we have presented and developed a 
hadron-level "meson emission model" of charmonium de- 
cays of the type '3/ — ?> NNm, where is a generic char- 
monium resonance, is a nucleon and to is a light meson. 
The model assumes that the decays take place through 
meson emission from the nucleon or antinucleon line, as a 
hadronic "final state radiation" correction to a 4' NN 
transition. As the model is relatively simple, we are able 
to evaluate the predicted DP event densities for many ex- 
perimentally accessible and measured processes; in par- 
ticular we give event densities for 5* = ?7c, J/ip (and ip'), 
XcOj Xci and ip' and to = tt'^, /□ and w, and implicitly all 
cases with the same J^^ quantum numbers. 

We used the reaction J/tp ^ ppn^ as a test case with 
no free parameters (the J/^pp and NNtt couplings are 
known), and compared the meson emission model pre- 
dictions for the projected event density in Mp^ and the 
partial width T{J /%p ppir^) to experiment; the partial 
width is in good agreement, and the Mp^ event density 
appears to describe the non-A^* "background" contribu- 
tion to this reaction observed experimentally. We also 
give predictions for r(^' — ppir'^) for all 5" cases consid- 
ered here; the overall trend of large and small widths and 
their approximate scale is reproduced by the model. 

We also considered scalar and vector meson produc- 
tion. We estimated 5* — > ppfo branching fractions for 
a light /o(500), and noted that the J/tp and tp' are fa- 
vored for these studies, and the ip' is favored for glueball 
searches. In vector production we considered J/ip — ppuj 
in particular, and noted that a high statistics study of 
this reaction could be used to estimate the NNuj cou- 
plings (g^j and K^j), which play a crucial role in meson 
exchange models of the NN force. We showed that the 
J/4' PP^ Dalitz plot event density is rather sensitive 
to the poorly known NNlo Pauli coupling k^. Determi- 
nation of meson-nucleon strong couplings is a potentially 
very interesting application of ^E* — > ppm decays. 

There are several theoretical developments that will 
be very important for future applications of this model. 
One should incorporate A^* resonances; this is conceptu- 
ally straightforward but may be technically complicated, 
as it will introduce many new and poorly known reso- 
nance coupling parameters and phases. This develop- 
ment is of course crucial to describe the data in reac- 
tions such as J /%p ppT^'^ , which clearly shows A^* res- 
onance peaks (Fig[3]). Another important development 
is the substitution of plausible 5r7V(*)iV(*) and N^*'^Nm 
hadron vertex form factors for the assumed coupling con- 
stants; the difficulty here is that hadronic form factors are 
poorly known, and models such as "^Pq that predict form 
factors have not been adequately developed and tested. 
Another interesting generalization of the strong vertices 
assumed here would be to include a J /ippp Pauli cou- 
pling; as we noted previously pi!], this can explain the ob- 
served e^e~ J /ip pp angular distribution. Finally, 
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one should include other significant decay mechanisms, 
as they become apparent through high-statistics studies 
of experimental Dalitz plots. These additional mecha- 
nisms might include intermediate meson resonances m! 
that couple strongly to NN , as in m'm — > ppm; 

if the m' resonances lie in the physical region, they will 
give rise to characteristic m' resonance bands in Mpp that 
could be identified and incorporated in a more complete 
decay model. 
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